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Abstract 

Geodesic equations of timelike and null charged particles in the Ernst metric are 
studied. We consider two distinct forms of the Ernst solution where the Maxwell 
potential represents either a uniform electric or magnetic field. Circular orbits in 
various configurations are considered, as well as their perturbations and stability. 
We find that the electric field strength must be below a certain charge-dependent 
critical value for these orbits to be stable. The case of the magnetic Ernst metric 
contains a limit which reduces to the Melvin magnetic universe. In this case the 
equations of motion are solved to reveal cycloidlike or trochoidlike motion, similar 
to those found by Frolov and Shoorn around black holes immersed in test magnetic 
fields. 
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1 Introduction 


The Ernst metric [1], also known as the Schwarzschild-Melvin metric or the electri¬ 
fied/magnetized Schwarzschild black hole, describes a black hole immersed in an uniform 
electric or magnetic field. In this paper we will distinguish the two cases by calling the 
former the electric Ernst spacetime and the latter magnetic Ernst spacetime. The solu¬ 
tion can be derived by applying a Harrison transformation [2] to a seed metric taken to 
be the Schwarzschild solution. The result is an exact solution to the Einstein-Maxwell 
equations parametrized by the black hole mass and the strength of the Maxwell field. The 
spacetime curvature is caused by both the black hole and the electromagnetic field. The 
Ernst solution was found just a year after Wald considered Kerr black holes immersed in 
a weak magnetic field [3]. In Wald’s solution a (test) magnetic field is found by solving 
Maxwell’s equation in a Kerr background. Hence the field in Wald’s construction is a test 
field which does not influence the spacetime curvature and can be considered as a limiting 
special case of the full Ernst solution. 

Despite being non-asymptotically flat, the magnetic Ernst solution is a useful model 
which describes black holes in certain astrophysical situations. Hence past work on 
geodesics focuses primarily on the magnetic case [4, 5]. It is likely that strong magnetic 
fields may exist in the vicinity of stellar mass and supermassive black holes. In particular, 
the motion of charged particles in this spacetime is studied in [6]. The geodesic motion 
for null and timelike particles is considered in [7, 8, 9]. Some later works on this topic 
include [10, 11, 12]. The motion of charged particles was also considered in the context 
of chaos and nonlinear dynamics in [13, 14], 

Nevertheless, in most realistic astrophysical situations the magnetic field surrounding 
a black hole is relatively weak such that the influence on the spacetime curvature is fairly 
negligible [15]. The motion of charged particles in this test field regime was studied in 
Refs. [16, 17, 18, 19], where the electromagnetic field does not influence the spacetime 
curvature and only affects motion of the charged particle via Lorentz interactions. Frolov 
and Shoorn [18], in particular, studied the case where the particle executes a curly cycloid¬ 
like trajectory which occurs due to the combination of a central gravitational force from 
the black hole on the particle together with a velocity-dependent outward Lorentz force. 
In this paper, we will also consider a similar motion in the Melvin spacetime, where, in 
the absence of a black hole, the inward gravitational force is provided by the magnetic 
field itself. 

The case of orbits around a black hole immersed in an electric field has received 
less attention compared to its magnetic counterpart. This case physically describes an 
electrically charged particle under the gravitational influence of the black hole and the 
axisymmetric electric field, in addition to a constant Coulomb force in the z direction. 
This can be viewed physically as a central force motion subject to an additional uniform 
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constant force, we should expect some similarities to particle motion in the vacuum C- 
metric [20]. 

This paper is organised as follows: In Sec. 2 we will review the Ernst metric and some 
of its relevant properties. In Sec. 3 we derive the equations of motion for a charged particle 
in the electric Ernst spacetime and study some solutions which are accessible analytically 
and perturbatively. Section 4 follows by deriving the equations of motion for the case of 
charged particle in the magnetic Ernst spacetime. For neutral particles the equations of 
motion for the electric and magnetic cases are identical. This includes the motion of null 
(lightlike) particles. This case has already been thoroughly studied in Refs. [7, 8, 9] using 
the Hamilton-Jacobi formalism. For completeness we consider them in Sec. 5 using the 
Lagrangian formalism we adopt for this paper. Some concluding remarks are discussed in 
Sec. 6. 


2 Ernst spacetime 

The four-dimensional Ernst metric [1] is given by 1 


ds 2 = A 2 (-fdt 2 + / Mr 2 + r 2 d9 2 ) + A 2 r 2 sin 2 9 d0 2 , 

f = 1-—, A = 1 + -B 2 r 2 sin 2 9, (1) 

r 4 


where m parametrizes the mass of the black hole and B parametrizes the strength of the 
uniform electric/magnetic held. 2 If the black hole is immersed in a magnetic held, the 
Maxwell potential is 


Am = 


r 2 sin 2 9 

B --— dtp. 

2A 


( 2 ) 


On the other hand, the Maxwell potential corresponding to an electrihed black hole is 


Ae = B fr cos 9 dt. (3) 

It is worth noting that with the presence of the electromagnetic held, the Ernst solution is 
not asymptotically hat. The Maxwell tensor for either case is given by F = d fJt A v ^d v A IJt . 
The metric (1), together with either (2) or (3) is a solution to the four-dimensional 
Einstein-Maxwell equations with zero cosmological constant. We can easily see that set¬ 
ting m — 0 reduces the solution to that of Melvin [22], while for B — 0 we recover the 

1 A review and related references can be found in [21]. 

2 We will keep the notation where B represents the field strength parameter for both the electric and 
magnetic Ernst solutions, as they will be discussed separately in this paper. Thus it will be clear from 
the context that B represents the electric field strength for the electric Ernst solution (Sec. 3) while it 
represents the magnetic field strength in the analysis of the magnetic Ernst solution (Sec. 4). 
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usual Schwarzschild solution. 

We also note that there are some useful symmetries of the metric. For the case of the 
magnetic Ernst spacetime, (1) and (2) are invariant under the transformation 

0 -> -0. (4) 

For the case of the electric Ernst spacetime, the solutions (1) and (3) are invariant under 

d^n-9. (5) 

With these symmetries, we can consider only the case B > 0 without loss of generality. 


3 Charged particles in the electric Ernst spacetime 


3.1 Equations of motion 

The motion of a test particle of charge per unit mass e is described by a trajectory x^{r) 
where r is an appropriate affine parametrization. In the case of timelike particles, r may 
be regarded as the proper time measured by the particle. The motion is determined by 
the Lagrangian C = ^g^x^x 1 ' + eA fl x tl 1 where overdots denote derivatives with respect to 
r. The subsequent equations of motion can be derived using the Euler-Lagrange equation 
57 Hi 1 ' = 77^ • F° r the case °f charged particles in the electric Ernst spacetime, where the 
Maxwell potential is given by Eq. (3), the corresponding Lagrangian is 


C 


1 

2 




y2 


+ 7 + r 


2n2 


r 2 sin 2 9 - 9 


+ eBfr cos 9t. 


( 6 ) 


Since d/dt and <9/<90 are Killing vectors, they give rise to constants of motion E and 
<f>, which we may interpret as the energy and angular momentum of the particle. The 
conserved quantities reduce the equations for t and 0 into first integrals: 


i = 


E + eBr cos 9 

A 


A 2< f> 
r 2 sin 2 9 


(7) 
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Applying the Euler-Lagrange equation to the remaining two coordinates gives 



( 8 ) 


(9) 


Here, the primes appearing in f indicate derivatives with respect to r. The invariance of 
g flu x ll x L/ = e, together with (7), gives a first integral equation 



( 10 ) 


By appropriately rescaling the affine parameter r, we can set the magnitude of e to be 
unity if it is nonzero. Therefore, for timelike particles we have e = —1, and e = 0 for null 
(massless) particles. 

We note that due to the symmetry of the solution shown in (5), we may assume 
without loss of generality that both B > 0 and e > 0. Since if both are negative the 
above equations remain unchanged, if they have opposite signs, the resulting equations 
are equivalent under (5). 

If we set m = 0 in the above equations, they reduce to the geodesic equations of 
the electric Melvin spacetime, while setting B = 0 reduces to the well-known geodesic 
equations aroimd a Schwarzschild black hole. Furthermore, setting B = b/e and neglecting 
terms of order 1/e 2 and beyond describes the motion around a black hole immersed in a 
test electric field. A simpler, though slightly less rigorous way to describe this is obtained 
by setting A = 1 in the above equations. In the test field regime, the electric field does not 
influence the spacetime curvature, hence its effect on the particle is purely electro dynamic 
in nature. 

3.2 Curves of zero velocity 

In the present case where the black hole is immersed in an electric field, it is not possible 
to cast Eq. (10) in the form of of an effective potential equation. Nevertheless, it is still 
possible to study the existence of bound and unbound orbits by rearranging to obtain 



( 11 ) 
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From the above equation, we can find regions accessible to the particle where the coordi¬ 
nates r and 6 lead to positive values in the right-hand side of Eq. (11). The boundaries 
of these regions are points where r = 6 = 0, and therefore represent turning points of 
the trajectory in which the particle reaches zero velocity. Hence, in earlier works such as 
Refs. [23, 24] call such boundaries curves of zero velocity. 

Figure 1 shows some typical examples of the regions accessible by the charged particle. 
The r and 6 coordinates are represented by the horizontal and vertical directions, respec¬ 
tively. The shaded areas represents the areas inaccessible to the particle. The top row is 
for the case $ = 3.5, while the bottom row are plots for <J> = 4. The leftmost diagrams 
on both rows show the case corresponding to B — 0, which is simply the Schwarzschild 
effective potential. As can be shown from the well-known Schwarzschild geodesics, around 
the value of $ = 4, the Schwarzschild potential contains a finite well, in which the particle 
is in a bound orbit around the black hole [25, 26]. When <f> is reduced to about to 3.5, 
the potential barrier disappears and the particle may fall into the horizon at r = 2m. 

The electric field is turned on and increased in the second, third and fourth plots 
of each row. As we see by visual inspection, the accessible regions are mostly pushed 
below 6 < 7t/ 2 (towards the north of the equator). This corresponds to the intuitive 
interpretation that the north-pointing electric field tends to push the positively charged 
particle along that direction. If the particle is negatively charged, or if the field points in 
the opposite direction, the particle will be pushed southwards, and the above diagrams 
should be reflected along the 6 = it/2 axis where the equations of motion are equivalent 
under (5). 


3.3 Nearly circular orbits around weakly electrified black holes 

When B is small, the spacetime can be regarded as a perturbation of the Schwarzschild 
solution. In this subsection we shall consider what happens to circular Schwarzschild 
orbits in the presence of a weak electric field. 

We first consider the case where the electric field is parallel to the orbital plane; hence, 
in our coordinate system this is a polar orbit, where = 0, and now we only have angular 
momentum in the 9 direction, which we denote as L. Recall that circular Schwarzschild 
orbits of radius ro have energy and angular momentum given by [25, 26] 


E 2 


ro — 2m 
r 0 (r 0 - 3 m) 1 


r^m 

r 0 — 3m 


The solution to the equations of motion is 


r = r 0 = constant, 9{r) 



( 12 ) 


(13) 
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Figure 1: (Color online) The accessible regions for charged particles in the electric Ernst 
spacetime of charge e = 1, E = 0.963. The vertical and horizontal axes correspond to 9 
and r, with the shaded regions as those inaccessible to the particle. The top and bottom 
row correspond to particles of angular momentum <f> = 3.5 and 4.0, respectively. In each 
row, from left to right the electric field strength is B — 0, 0.0005, 0.0010, and 0.0015. 
The above diagrams are plotted for the range r e [2m, 27m] and 6 e [0,7r]. The equator 
6 = 7t/ 2 is represented in each diagram with the horizontal dotted line. 

For small B , suppose the trajectory of the particle is given by 

r{r) = r 0 + Bri(r) + O (B 2 ) , 0{r) = + B9 1 (t) + O (B 2 ) . (14) 

r o 

Substituting into Eqs. (8) and (9) and solving the equations to linear order in B , we 
obtain 3 


where 


r i 


- u> 2 ri + 


e(3r 0 — 8m) 
V r o( r o - 3m) 


cos h2r, 


2 _ m(r 0 - 6m) 2 

^ r o( r o — 3m) ’ 


2 

m 

r o ( r o — 3m) 


(15) 


(16) 


Similar to the case of orbits around weakly accelerated black holes [20], the radial equation 
reduces to that of a harmonic oscillator with a periodic driving force. Stable orbits 
corresponding to cu 2 > 0 are satisfied by r 0 > 6m, in accordance with the well-known 
results regarding the stability of Schwarzschild orbits. For the case of charged particles 
with e ^ 0, there is an extra term which acts as a periodic driving force. 

3 The procedure is similar to the perturbations considered in [20], where the details can be found. 
Additionally, it is important to note that we have the implicit assumption that e is sufficiently small that 
terms of the order O ( eB 2 ) may be neglected. 
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We can check that the full nonperturbative numerical solutions of (8) and (9) are 
consistent with the above results if they are solved for small B. Figure 2a shows the 
numerical solution of r vs r of a neutral particle in nearly circular polar orbit around the 
black hole. It is possible to check that the particle undergoes simple, undriven oscillation 
in the radial direction with a period of 27r/c o ~ 233, in agreement with Eq. (16). 

Figure 2b shows a similar solution for a charged particle of e = 1, demonstrating the 
behavior of a driven oscillator with driving frequency given by 12. 




(a) r vs t for e = 0. (b) r vs r for e = 1. 


Figure 2: (Color online) Plot of r vs r of perturbations about circular polar orbits of 
r 0 = 7 (in units where m = 1) by a weak electric held B = 0.0001. The plot shows the 
numerical solutions of r vs r for (a) a neutral particle e = 0, and (b) particle of charge 
e = 1. From Eq. (16) we can see that the period of oscillation in (a) is 2tt/uj ~ 233, in 
agreement with the oscillations seen above. In (b) there is an additional driving force 
with frequency 12. The initial conditions are r(0) = 7.05, r(0) = 0, and 0(0) = 0. 


Next we we consider the case where the initial circular orbit is in the equatorial plane. 
In this case the perturbed equations decoupled and are 


r i = — cu 2 ri, 


9l v/r„(r„-3 m)’ 


(17) 


where oj and 12 are the same as defined in (16). This time we see that the equation for 
61 has a constant term if e ^ 0. This represents the constant force by the electric field, 
pushing the particle out of the equatorial plane. 


3.4 Circular orbits in arbitrary field strengths 

For B that is not necessarily small, we can End circular orbits by demanding that r and 
6 be constant in Eqs. (8) and (9). This requires that the energy and angular momentum 






be given by 


E 

cD 2 


eBro (ro — 2 m — m cos 2 6*o) 
m cos 6*o 

256e 2 £> 2 rQ sin 4 6 0 (4m + 2B 2 r^ — 3mB 2 r q — mB 2 rl cos 2 0 O ) 
m 2 cos 2 (4 — B 2 rl sin 2 0 O ) 5 


for some constant r 0 and 0 O - 

We perturb about the circular orbits by writing 


( 18 ) 


r(r) = r„ + er^r) + O (e 2 ) , 9(t) = 6„ + eO^t) + O (e 2 ) . (19) 

Substituting into Eqs. (8) and (9) and expanding to first order gives 



where A t] are somewhat complicated functions of r 0 and 9 0 , though it can be handled 
appropriately with the aid of a symbolic computation software such as MAPLE. We find 
the normal mode frequencies by solving for the eigenvalues of the above 2x2 matrix, 
which is given by 

A± = — [a u + A22 ± y /(An — A22) 2 + 4^21^12^ . ( 21 ) 

The orbit is unstable if one or both of the eigenvalues are positive. Hence, a sufficient 
condition for instability is that the larger eigenvalue has a range of ro which is positive. 
The left-hand plot in Fig. 3 shows A + as a function of ro for fixed e = 1. When the held 
strength is relatively weak, there is a wide range of r 0 where stable orbits are possible. 
This range shrinks as mB is increased. Beyond the value of mB clit ~ 0.00392360, there 
are no more stable orbits. 

The particle’s charge is another free parameter and hence /i cn) should also depend on 
e. The right-hand plot of Fig. 3 shows the dependence of mB clit on the charge e. We can 
clearly see that the greater the particle’s charge, the lower the value of mB cvlt . Indeed, 
when the particle’s charge is higher, it will experience a greater Coulomb force from the 
electric held and hence, be more easily rendered unstable. 
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Figure 3: (Color online) Left: Plot of A + vs r 0 , of various mB, with fixed e = 1 and 
m — 1. From the bottom curve, the values of B are 0.001, 0.003, 0.0039236(= £> crit ), and 
0.001. Right: Plot of B ciit vs e. 


4 Charged particles in the magnetic Ernst spacetime 

4.1 Equations of motion 

For the case of the magnetized Ernst spacetime, the vector potential is given by (2) and 
the Lagrangian for the charged test particle is 


c = \ 


A 2 (-/f 2 + y+r 2 0 2 


r 2 sin 2 0 ■, 


+ eB 


r 2 sin 2 6 • 
2A 


( 22 ) 


In this case, d/dt and d/d(j) are still Killing vectors associated with the conserved energy 
and angular momentum E and $. This leads to the first integrals, 


t = 


A 2 


A 2 /’ r 2 sin 2 # 

The equations of motion for r and 6 are 


$ _ 


eBr 2 sin" 9 
2A 


(23) 


r = 


r 2 + / I r + 


f__dA 

2/ A 

/ 

r 3 sin 2 6 


eBf ( rd r A 
1 ~ ~2A~ 


r 2 d r A 


• 2 2d e A^ fd r A , f \ E 2 


6 Z - 


1 - 


rd r A 


A 


e = 


rA 
dftA ( r 2 


$ 


A ) A 

eBr 2 sin 2 6 
4 -— 

eBr 2 sin 2 6 
A 


-rd — 


A + 2/7 A 4 


(24) 


A 1 

Wf~ 

1 

1 ( 

1 r 4 

sin 3 9 { 

+ 

eB .i 


- e 2 - 2 i + 


d r A\ rd E 2 doA 


r 2 A sin 9 


cos 9 — 
cos 9 - 


sin 9DqA 

A 

sin 9doA 
2A 


A 1 r r 2 A 5 f 

2 -„2 ,/ i \ 2 


$ _ 

$ - 


eBr 2 sin" 9 
2A 

eBr 2 sin 2 6 1 
2A 


(25) 
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and the first integral now is 


E 2 9 1 ( eBr 2 sin 2 9 

-h A - $- 

A 2 / r 2 sin 2 # \ 2A 


+ J +r 


2n2 


= e. 


(26) 


As expected, setting B — 0 reduces to the geodesic equations of the Schwarzschild space- 
time, while m — 0 gives the equations for the Melvin magnetic universe. Setting A = 1 
describes charged particles in the test field regime considered in [18], where the mag¬ 
netic field is sufficiently weak so as not to influence the spacetime curvature but only 
the charged particle through Lorentz interaction. A common special case of interest is 
equatorial motion where 6 = n/2 and 6 = 0. Then we have 6 = 0 and the motion is 
confined to the equatorial plane. Thus, the r equation can be solved by direct integration 
of Eq. (26). Furthermore, due to the symmetry given in (4), by a reasoning analogous to 
the electric Ernst case, we shall assume B > 0 and e > 0 without loss of generality. 

Equation. (26) can be expressed as an effective potential formulation in the following 
form, 


A 4 


r 2 


+ r 2 f6 2 ) = E 2 — V 2 S , V? s = A 4 / 


r 2 sin 2 6 


$ _ 


eBr 2 sin 2 6 
2A 


e 

A 2 


(27) 


4.2 Circular orbits in the test field approximation 

In the context of astrophysics, it is typically sufficient to consider magnetic fields which 
are sufficiently weak so that they do not influence the spacetime curvature [15]. In our 
present paper we shall call this case the “test field regime”, where, lacking the influence 
from the magnetic field, the spacetime metric is simply the Schwarzschild metric [A —^ 1 
in Eq. (1)], but the vector potential associated with the magnetic field is still given by 
( 2 ). 

Previous works such as [17, 16, 18] considered the motion of charged particles in the 
test field regime. The only influence of the magnetic field on the charged particle is via 
the Lorentz interaction. This regime is contained within our present equations of motion. 
The resulting equations of motion that follow from the test field regime can be easily 
obtained by setting A —> 1 in Eqs. (23)-(26). However we wish to have a more rigorous 
description by seeking an appropriate parameter that will allow a more explicit transition 
from the test field regime from the full Ernst metric. 

We begin by making the following observation: based on the order-of-magnitude es¬ 
timations of [18], even in the test field regime a charged particle may still experience a 
significant Lorentz force that influences its orbit. In other words, B is sufficiently small 
such that it does not produce spacetime curvature (A ~ 1), but the Lorentz interaction 
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is still present; hence, eB is non-negligible. With these considerations, we write 

B = K (28) 

e 

and expand Eqs. (23), (24), (25) and (26) in powers of 1/e. At zeroth order, we have 
B — y 0, leading to A —» 1, but the Lorentz interaction eB = b is still present. In this 
order we recover the test field regime considered in [18]. Terms of the order 0(1/e) and 
higher come in as the gravitational correction due to the spacetime curvature induced by 
the magnetic field. 

As an example, we consider circular equatorial orbits in the magnetic Ernst spacetime. 
Taking r = r 0 = constant and 6 = 7t/ 2 in the full equations of motion (24)-(26), we find 
that the energy corresponding to a circular orbit is 

^2 ( r o — 2m) 2 (4e 2 + rib 2 ) 3 (44>e 2 + Qr^b 2 — 2e 2 brl) (^r^b 2 — 4<f>e 2 — 2e 2 brl) 

256e 8 rQ (3rlmb 2 — 4me 2 — 2rQ& 2 ) 

The corresponding values of $ may be found by substituting (29) into (26) and solving 
the quadratic equation for $, for which there are two solutions 


where 


= 


2e 2 r 0 


ro(3ro — 5m) b 3 + 4 e 2 r 0 mb ± 2 \[K 


(4e 2 + rib 2 ) [ro(3r 0 — 5 m)b 2 + 4e 2 (r 0 — 3m)] 


(30) 


K = 4rle 2 [(r 0 — 2m) 2 e 2 + 2 rl — 12r 0 m + 14m 2 ] b 2 — rp(3r 0 — 5m)(2r 0 — 3 m)b 4 

+ 16e 4 m(r 0 — 3m). (31) 

For concreteness, we analyze the innermost stable circular orbits (ISCOs) in the test 
field regime. These are orbits which are marginally stable, where all circular orbits with 
radius less than the ISCO are unstable and those with radii larger than it are stable. 
The stability of the circular orbits is determined by perturbing about the circular orbit 
solution by writing 


t (t) = r o + er 1 (r). 

With this ansatz, perturbing Eq. (24) to first order in e reduces to 

'f\ — — A r 1 , A = kiE 2 + /c 2 $ 2 + + k Al 


(32) 


(33) 
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Figure 4: (Color online) Plots of ISCO radii vs $ = <f>±, in units where m = 1. The lower 
(solid) curve represents the ISCO radius for $ = <f> + , while the upper (dotted) curve 
corresponds to $ = <f>_. 


where 

512 [(31mr 0 — 9 Tq — 27m * 2 ) r^b 2 + 4 (tq + 14m 2 — 10mr 0 ) e 2 rg6 2 — 16me 4 (r 0 — m)] 

1 rg(4e 2 + rg6 2 ) 6 (r 0 — 2m) 2 

16e 2 rg6 2 (ro — 2m) — rg(3r 0 — 8m)b 4 + 16e 4 (3r 0 — 8m) 

k2 = (4e 2 + r 2 b 2 ) 2 

4 b 3 [(3rg — 8 r 0 m)b 2 + 8e 2 ] e 2 

3 = (4e 2 + r 2 6 2 ) 3 ’ 

4 b 2 [r 0 5 2 (8m — 3r 0 ) + 4e 2 ] 

4 = (4e 2 + r 2 b 2 f ’ 

and E and $ are given by (29) and (30). As usual, stable circular orbits have A > 0 and 
unstable ones have A < 0. The expression for A is expanded in powers of 1/e to give 

A± ~ Af + + O (1/e 4 ) , (35) 

where the ± signs correspond to the two distinct choices of $ = To find an ISCO, we 
solve A± = 0, where the upper and lower signs lead to two distinct ISCO radii r 0 = r±. 
At zeroth order in 1/e, which is the test held regime, the results of [18] are reproduced, 
shown in Fig. 4. Particularly for b ^ 0, we typically have r± < 6m, i.e., the typical ISCO 
radius for nonzero magnetic fields is smaller than the Schwarzschild ISCO radius. 

Next we introduce the gravitational corrections to the above results. For ease of 
exposition, we shall denote the ISCO radius evaluated in the test held regime as the 
“test held ISCO radius”, and the ISCO radius when taking hrst-order corrections into 
account as “gravitationally corrected ISCO radius”. For an ISCO with given r± and b, 
we find the gravitational correction to A by calculating the second term in (35). For the 
background ISCO parameters, by their definition of being marginally stable orbits, we 
have A±^ = 0. Thus, if Ai 2) is found to be positive, the ISCO orbits are stabilized by 
gravitational corrections. Conversely, if we find A[)_ to be negative, then the ISCO orbits 


(34) 
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(a) A±^ vs r±. 



(b) A±^ vs b. 


Figure 5: (Color online) Plots of A® against r and b, plotted in units where m — 1. The 
solid curves represent A + and r + while the dotted curves correspond to A_ and r_. 


( 2 ,') 

are rendered unstable. The particular expression of A± is again too complicated to be 
displayed here, but the dependence of b or, equivalently, r± is shown in Fig. 5. 

We note that in the case <3>_, we have X_ >0; hence, circular orbits in the test held 
ISCO radius are stabilised by gravitational corrections. In other words, test held ISCO 
radius r_ is now stable and the gravitationally corrected ISCO radius is slightly smaller 
than the test held ISCO radius. 

The case <h + is more interesting. As seen in Fig. 5b, for relatively small b , we still have 

( 2 ) 

A+ > 0; the gravitationally corrected ISCOs are also further inwards. However, beyond 

( 2 ) 

a critical value of approximately b ~ 0.22134, the value of A+ is negative. Hence, in this 
case gravitational correction actually destabilizes the test held ISCOs. When gravity is 
taken into account, the gravitationally corrected ISCO radius is further away from the 
black hole. 

4.3 Curly orbits: Trochoidlike trajectories in the Melvin space- 
time 

In Ref. [18], Frolov and Shoorn considered an interesting behavior where the trajectory 
of charged timelike particle curls up into a cycloidlike (or more generally, trochoidlike) 
motion. This behavior is due to the fact that, for motion in the equatorial plane, the 
charged particle experiences a gravitational force in the direction orthogonal to the mag¬ 
netic held lines. As explained in the previous subsection, they have considered the test 
held (A —> 1) regime, where only the Schwarzschild black hole contributes to the or¬ 
thogonal gravitational force, and the magnetic held is solely responsible for the Lorentz 
force. 

Here we consider a similar effect for the Melvin magnetic universe. 4 While there is 

4 We expect that the presence of a black hole does not add any notable physical difference to bound 
trochoidlike orbits, aside from a stronger gravitational force towards the center. Hence we avoid unnec¬ 
essarily cumbersome equations by taking m = 0. 
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no black hole present in this case, the magnetic field itself exerts a gravitational force on 
the particle, in addition to providing the Lorentz interaction. The curling-up behavior is 
characterized by the fact that (j) changes sign in Eq. (23). This sign change occurs when 
r = r*, where 

1 

r* =- 

sin 9 

In the following we will consider only bounded motion in which the particle is confined 
within the range r min < r < r max , where the boundaries are defined by E 2 = V 2 E in 
Eq. (27). If we further consider motion confined in the equatorial plane, it is possible 
derive approximate solutions representing the cycloid- or trochoidlike motion if we consider 
perturbations about circular orbits in the Melvin spacetime. Hence, we solve Eq. (24) by 
substituting 9 = 7t/ 2 and taking 

r(r) = r 0 + eri(r) + O ( e 2 ) . (37) 


4$ 


2 eB - $H 2 ' 


(36) 


At zeroth order, demanding that r 0 is constant gives 

2 _ (16$ 2 - 4e 2 £> 2 rQ + AeB 3 $r* - B A $ 2 r*) (4 + B 2 r 2 ) 

~ 512r^B 2 ' 

At first order, the equations of motion reduce to fq = — cu 2 7q, where 

2 _ 2 (3$ 2 H 6 r§ - 12eB 5 $r% + 12 e 2 B 6 r$ - 16H 2 $ 2 r 2 + 128$ 2 ) 
U (4 + B 2 r 2 ) 3 rQ 

Substituting (37) into the 0 equation, we find 

(j) = a 0 + aqe cos cur + O (e 2 ) , 


(38) 


(39) 


(40) 


where we have used the solution r\ = cos cur, and 


(4$ + B 2 $r 2 - 2 eBr 2 ) (4 + B 2 r 2 ) _ B^r* - 2eB 3 r% - 16$ 

a 0 — --j a i — - 


167-q 


8 rg 


(41) 


Therefore, the approximate solutions to the equations of motion are, to first order in e, 


r(r) = r 0 + s cos cur, 
c p(t) = aoT -|-sm cur. 

UJ 


(42) 


This describes the locus of a trochoid, where there are three possible types depending on 
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the relationship among its parameters. Let us define 


OL\E 

If rj — 1, the motion is that of a usual cycloid, where the trajectory forms sharp cusps 
at maximum r. In more general cases, rj < 1 is known as the prolate cycloid, and 77 > 1 
corresponds to a curtate cycloid. 

We can verify the above solutions with the numerical solutions of the fully 11011 - 
perturbative equations with the appropriate range of parameters. Fig. 6 shows a solution 
for B = 0.05, $ = 5 and e = 16. It appears that, nearly circular orbits with r = r* lying 
close to the vicinity of the oscillation typically occur for highly charged particles; hence, 
the choice e = 16 in Fig. 6 . 

Of course, the above calculations are valid for nearly circular orbits where the radius 
is close to ro- Moving beyond the perturbative range, we can obtain the curly orbits from 
numerical solutions. Some examples are shown in Fig. 7. 

For the rest of this section we briefly consider the critical, cycloidlike orbits in the 
nonperturbative range. This class can be easily studied since they can be characterized 
by simple equations. It is clear that the trajectories develop cycloidlike cusps if r* coincides 
with r max . Since by definition V e g (r = r max ) = E 2 , it follows that 


( 43 ) 


r* 


2e(E-l) 

B ’ “ 'BE 


(44) 


We conclude that such orbits exist under the condition E > 1, and > 0. (Recall that 
due to Eq. (4), we may assume B > 0 and e > 0 without loss of generality.) Following 
the spirit of [27], one may further understand the structure of the orbits by identifying 
the periodic orbits since any generic orbit appears like perturbations of periodic ones. 
Periodic orbits are defined as trajectories which return precisely to their initial conditions 
after a finite proper time; thus, the orbits periodically retrace the same path repeatedly. 

Since, for equatorial orbits, we have 9 = n/2 and 6 = 0, we can take Eq. (27) as the 
first integral for the r motion. Together with the equation for 0 from Eq. (23), we can 
obtain a trajectory analytically by integrating 


r = cb = W£-;W-fr ) 2 + ;g .... 

0 d 0 A 2 ($-^) 

We can classify the periodic orbits by the number of cusps n. formed before the particle 
returns to the initial conditions. A particle moves between cusps by starting from r max , 
reaching a turning point at r m ; n , then returning to r max again. The orbit will be periodic 
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3.56 



3 52 - 1 - 1 - 1 - 

- 0.1 - 0.05 0 0.05 0.1 

X 


(a) E = 1.0079, ry ~ 0.608 < 1. 



(b) E = 1.007874016, 7/ ~ 0.997 ~ 1. 



x 

(c) E = 1.00785912, 7 / ~ 5.329 > 1. 

Figure 6 : (Color online) Numerical solutions for B = 0.05, $ = 5, e = 16. The dashed 
lines are the boundary of the allowed motion defined by E 2 = V 2 ^ while the dotted curve 
represents r = r*. In (a) we show the case corresponding to a prolate cycloid with 77 < 1, 
(b) is the common cycloid, and (c) is the curtate cycloid. As we see in (b) the curve 
r = r* coincides with the maximum allowed r of the motion, producing the cycloidlike 
trajectory. For each solution, for a given set of parameters the value r 0 is 

calculated from d (V^ s ) /dr = 0 and £ may be calculated from E 2 = V 2 S . 


17 




























(a) (B, e, E, $) = (0.05, 5, VL07,5) (b) (B, e, E, $) = (0.1,1,1.78,5) 




(c) ( B , e, E, $) = (0.2,1,1.6,5) (d) (B, e, E, $) = (0.05,1,1.2,5) 

Figure 7: (Color online) Trochoidlike trajectories for various orbital parameters. The 
dashed circles indicate the boundary of the allowed motion defined by E 2 = V 2 S while the 
dotted circle represents r = r*, the turning point of </>. In (c) it can be seen that r* lies 
beyond the accessible range and, therefore, does not curl up. 
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(a) n = 1, E = 4.6151 


(b) n = 2, A = 1.69857 




(c) n = 3, E = 1.2767 


(d) n = 4, £ = 1.0882 


Figure 8: (Color online) Periodic cuspy orbits for B = 0.3, and e = 1. 


if the difference in 0 between the cusps is some rational fraction of 2n. In terms of (45), 



for some integer n. In practice, one can find a periodic orbit in the following way: Given 
parameters ( B,e,E ), the angular momentum <h and r* = r max are determined from (44), 
while r m i n is determined from solving E 2 = V 2 S . With these parameters one can find A0 
from Eq. (46). Periodic orbits are found by tuning one of the parameters, say, E for some 
fixed B and e to find an integer solution for which A0 = 2n/n. In Fig. 8, we show an 
example of periodic orbits for n ranging from 1 to 4. It appears that for fixed B and e, 
the number of cusps increase as the energy E decreases. 
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5 Neutral particles 

5.1 Effective potential for neutral particles 

For the case of neutral particles, the trajectories are then purely determined by the space- 
time curvature, and they are not affected by the electromagnetic held. Hence for e = 0, we 
see that the geodesic equations for electric and magnetic Ernst spacetime becomes iden¬ 
tical. In particular, Eq. (10), now identical to (26) can be cast into an effective potential 
form 

A* (*’ + = £2 - ^ = A V - £) ■ (47) 

To see the structure of the effective potential for photon orbits, we set e = 0 in Eq. (47) 
and find curves in the (r,6) plane which satisfies E 2 = V 2 S . These curves correspond to 
points where r = 9 = 0 and serve as the boundary of the region accessible by the photon. 

Figure 9 shows a typical structure of the effective potential of photon orbits in the 
Ernst spacetime, where the shaded region indicates areas not accessible to the photon of 
a given E and <3>. The first diagram from the left is the plot for B = 0, which is simply the 
well-known effective potential for Schwarzschild orbits. The following diagrams towards 
the right show the effect of increasing B for fixed E and <3>, where we see the “neck” 
gradually pinches off as a potential barrier forms, creating an isolated finite potential well 
in the third diagram. The potential well vanishes if B is increased further. The critical 
value of B where the potential value shrinks to a single point corresponds to the case of 
circular photon orbits of constant r and 9 , which we will study in detail in the following 
subsection. 

Figure 10 shows the effective potential for neutral timelike particles. The effect is less 
interesting in this case. The first diagram on the left is the potential for Schwarzschild 
orbits with B = 0. Turning on the field strength binds the orbiting particle more closely 
to the center, hence shrinking the accessible region of the particle. 

5.2 Stability of circular photon orbits 

It is possible to solve the equations of motion to End circular photon orbits [7] where r 
and 9 are constant. To hnd such solutions we begin by first observing that for e = 0, 
there is a simple solution to the 9 equation with 9 = n/2. It follows that substituting 
9 — 7r/2 in (47) reduces it to a one-dimensional radial potential equation, where circular 
orbits can be found by considering a value of r = ro which satisfies 

T[V&(r,*/2)]=0. (48) 


20 




Figure 9: Effective potential for photon orbits for m — 1, E = 0.669 and $ = 3. From 
left to right, the field strength is B = 0, 0.17, 0.175 and 0.18. The shaded regions indicate 
areas where V 2 S > E 2 and, hence, not accessible to the particles. Horizontal and vertical 
axes respectively, correspond to r and 9. The horizontal dotted line through the center 
indicates the equator at 9 = tt/2. The ranges for the above plots are r E [2,10] and 
9 E [0,7r]. 





Figure 10: Effective potential for neutral timelike orbits for m — 1, E — 0.97 and $ = 4. 
From left to right, the field strength is B = 0, 0.01, and 0.03. The shaded regions indicate 
areas where V^ s > E 2 and, hence not accessible to the particles. Horizontal and vertical 
axes respectively correspond to r and 9. The horizontal dotted line through the center 
indicates the equator at 9 = 7t/ 2. The ranges for the above plots are r E [2,20] and 
9 E [0,7r]. 
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Thus it follows from (48) and (47) that the equation for circular, equatorial orbits of 
radius r = ro requires the field strength and energy to be given by 


l2 = 4(r 0 - 3m) = /(r 0 )A(r 0 , f ) 2 ^ 2 

r o(3r 0 —5m)’ Tq 


(49) 


It appears that $ is a free parameter, provided that E and B satisfies (49). As shown 
in [28], for a given value of mB , there are two possible radii in the range r$ > 2m which 
satisfies (49). The inner radius is unstable while the outer one is stable. The inner and 
outer radii coalesce at the critical value of 


ui B 2/3 crit , 


(50) 


where 


ft crit — 


3 (a/19 — 1) 


\ (\/l9 + 8) 2 (\/l9 + 3) 


~ 0.094 683 19. 


(51) 


For values satisfying mB > /3 crit , there are no circular orbits. For a given m and B 
where mB < 2/3 cr it, there are two circular orbits, where the one with the smaller radius is 
unstable and the outer one is stable. When mB is tuned to be equal to 2/3 cr i t , the inner 
and outer orbits coalesce and the circular orbit is marginally stable. The numerical value 
of /3 CI ^ was first obtained in [7] by numerical root-finding, while the analytical expression 
was obtained later in [28] by inspecting the condition given in Eq. (49). 5 We will now 
demonstrate the stability of the orbits explicitly in a calculation which also takes into 
account possible motion in the 9 direction. Perturbing around the circular orbit solution, 
we write 


r(r) = r 0 + £r 1 (r), 9 (t) = | + e0i(t) + O (e 2 ) , (52) 

where £ is a small perturbation parameter. Substituting (52) into Eqs. (8) and (9) [or 
(24) and (25)] and expanding, we find that the equations for r and 9 decouple at linear 
order, 



5 To compare the results, it is worth noting that the definition of B in [7] differs from (1) by a factor 
of 2. 
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(a) r = ro — O.OOl, ro = 10. 


(b) r = ro — 0.001, ro = 4. 


Figure 11: Plots of r vs r of perturbed circular photon orbits about (a) r 0 = 10 and 
(b) r 0 = 4; in both cases, the angular momentum is $ = 2. The equations of motion 
are solved in units where m — 1; for case (a) we have u; 2 > 0; hence, it shows a stable 
oscillation about r 0 = 10. We can check that the numerical solution agrees with the 
analytical approximation T = 2ti/u j r ~ 319. In case (b), r 0 = 4 leads to u; 2 < 0, the 
perturbed orbit falls into the horizon. 


where 


OL = 


3r 0 2 


16mr 0 + 15m 2 


2rg(r 0 — 2m) 




m 


UJn = 


r o( r o — 2m) 




(54) 


The stability of the orbits is ensured if both normal mode frequencies ay and ug are real. 
As we can see, Mq is always positive while in the r direction, while a; 2 is negative for the 
range 


y/l9 


m < 2m < r 0 < 


+ \/l9 


-m. 


(55) 


where the orbits are unstable. Therefore, circular photon orbits are stable for tq > 
8+ ^ m. The critical value of ro = 8+ ^ m corresponds with having mB = /3 cr it in 
Eq. (50), the value where the inner and outer radii coalesce. Thus, we have demonstrated 
explicitly that the inner radii are unstable. Furthermore, we have the additional result 
that the orbits are always stable in the 6 direction. 

The calculations can be verified by the numerical solutions. For example, perturbing 
about a circular photon orbit of radius r 0 = 10m results in an oscillation with a period of 
T ~ 319, shown in Fig. 11a in units where m = 1. This agrees with T = 2n/(jO r where ay 
is given by (54). Fig. lib shows an unstable orbit in the range given by 55, specifically 
r 0 = 4m. The instability of such an orbit is clearly seen as the perturbed particle falls 
beyond the black hole horizon at ro = 2m. 

Moving beyond perturbed circular orbits, we End by numerical integration more gen¬ 
eral bound orbits with nonconstant r and 6. By adjusting 5, E and $, we can hnd specific 
parameters where there exists a finite potential well. One such case is shown in the third 
figure in Fig. 9. An example of an orbit is shown in Fig. 12. 
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Figure 12: Photon orbit (e = e = 0) bound in a finite potential well. The parameters of 
this solution are B = 0.105, E = 0.6, and $ = 4, given in units where m — 1. 


5.3 Stability of circular timelike orbits 

The perturbation equations for a timelike circular orbit are more complicated than for the 
photon case. Therefore, we do not have a simple closed-form expression for the stability 
of orbits analogous to (51), and our result still requires numerical root-finding already 
considered in [8]. 

Thus, we will only consider this case briefly for completeness. For a given and B , 
circular timelike orbits have energy and angular momentum given by 


$ 2 = 


16?q ( 2B 2 r o — 3.B 2 m?Q — 4m) 


E 2 = 


(4 + B 2 r 2 ) (5 B 2 r 2 m - 3 B 2 r% - 
(r 0 — 2m) (4 + B 2 Tq) ( 2B 2 rlm 


12m + 4 r 0 ) 


B 2 rl 


8m + 4 r 0 ) 


r 0 (5B 2 r^m — 3B 2 Tq — 12m + 4r 0 ) 


(56) 


Perturbing around the orbits in the case of timelike orbits also yields decoupled equations 
of the form (53), but with the frequencies given by 


= 


16 


[(127-q — 37mr 7 0 + 30m 2 ro)5 6 


rj] (4 + B 2 Tq) (5£> 2 mrQ — 3 B 2 Tq — 12m + 4r 0 ) 


= 


(204mrg — 200m 2 rQ — 48 r^B 4 + (672m 2 rg + 128rQ 
624rQm)5 2 + 64mr 0 — 384m 2 ] 

16m (4 — 5 2 Tq) 


Tq (4 + B 2 Tq) (5B 2 rl — 3 B 2 r^ — 12m + 4r 0 ) 


(57) 

(58) 


Stable orbits require cu 2 > 0 and ujg > 0; the value of r 0 required for stability differs for 
different B, unlike in the photon case. 
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6 Conclusion 


In this paper we have considered the geodesic equations for charged particles in the Ernst 
metric. The metric represents a Schwarzschild black hole immersed in an axisymmetric 
electric or magnetic field. In the case of the electric Ernst metric, the geodesic equations 
describe charged particles experiencing a central force in addition to a uniform constant 
force along the axial direction. We find that nearly circular polar and equatorial orbits 
of radius r 0 > 6 m are stable in weak electric fields, in accordance with the stability of 
Schwarzschild circular orbits. In the case of polar orbits, the radial equation reduces 
to that of a harmonic oscillator driven by a periodic driving force. This driving force 
corresponds to the periodic frequency of the motion in the 6 direction. When the held 
strength is not necessarily small, it is still possible to have circular trajectories whose 
orbital plane is parallel to the equatorial plane but lies at a fixed distance above it. 

For the case where the Maxwell held is purely magnetic, we considered curly trajecto¬ 
ries in the Melvin spacetime. Such orbits may form if the charged particle simultaneously 
experiences an inward gravitational force which is counteracted by a Lorentz force di¬ 
rected outwards. Such motion exists already in the domain of classical electromagnetism 
for particles in uniform crossed electric and magnetic helds [29]. For particles in the 
Melvin spacetime, instead of an electric Coulomb force we have the geodesic motion in a 
spacetime curved magnetic held itself. 

We have also considered orbits in weak magnetic helds by a perturbative expansion of 
the full equations of motion in the Ernst spacetime. With this expansion, we recover the 
equations of motion in the Wald’s construction of weakly magnetized black holes, where 
the magnetic held is a test held that does not influence the spacetime curvature. By taking 
higher-order terms of the expansion, we are able to calculate gravitational corrections to 
the results of the test held case. 

Neutral particles do not distinguish between the magnetic and electric nature of the 
helds; hence, the equations of motion for e = 0 for the electric and magnetic cases reduce 
into each other identically. This case is studied in [11] and [8] by focusing on motion 
confined on the equatorial plane. Here we have shown explicitly that photon orbits of 
radius r 0 > 8+ ^ m are indeed stable even when the particles are perturbed slightly away 
from the equator. 

In this paper we have focused exclusively on electrically charged and neutral timelike 
particles, in addition to (neutral) photons. We have not explicitly considered more exotic 
particles such as magnetic monopoles or charged massless particles. Theoretically, the 
motion of magnetic monopoles in the magnetic Ernst spacetime should be identical to 
electric monopoles in the electric Ernst spacetime. So we should expect most of the results 
obtained already for electrically charged particles should carry over after performing the 
appropriate duality operations. The motion for charged massless particles is contained in 


25 



the equations of motion in Sec. 3 and 4, by considering e ^ 0 and e = 0. 
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